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Abstract: We report on the experimental and theoretical study of the 
spatial fluctuations of the local density of states (EM-LDOS) and of the 
fluorescence intensity in the near-field of a gold nanoantenna. EM-LDOS, 
fluorescence intensity and topography maps are acquired simultaneously 
by scanning a fluorescent nanosource grafted on the tip of an atomic force 
microscope at the surface of the sample. The results are in good quantitative 
agreement with numerical simulations. This work paves the way for a full 
near-field characterization of an optical nanoantenna. 
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1. Introduction 

Much interest has recently been paid to plasmonic nanostructures due to the capacity they of- 
fer to enhance light-matter interaction of elementary dipoles such as fluorescent molecules and 
quantum dots. This encompasses three different mechanisms which are often hard to disen- 
tangle. (1) Light absorption can be enhanced by an optical antenna, leading to an increased 
effective absorption cross section. This can be advantageously used in photodetection and pho- 
to voltaics JT|. (2) Plasmonic nanostructures may also induce significant changes in the sponta- 
neous emission dynamics by the Purcell effect. As an example, large Purcell factors associated 
to electromagnetic modes strongly localized in subwavelength volumes have been reported at 



the surface of disordered thin metallic films, which opens new perspectives such as the inves- 
tigation of strong coupling in a regime where high Purcell factors coexist with high absorption 
El El El- (3) A change of fluorescence intensity is also expected when a dipole is in near- field 
interaction with a plasmonic nanostructure. In an optical antenna, a field enhancement occurs, 
but it is often obscured by non-radiative processes leading to fluorescence quenching. 

The local density of electromagnetic states (EM-LDOS) is the basic quantity which governs 
these three mechanisms. As it rules the energy stored in all available modes, it provides a direct 
measurement of the probability of spontaneous light emission. The decay rate of a fluorescent 
dipole is proportional to the EM-LDOS. While in vacuum the EM-LDOS is homogeneous, it is 
known that it can be significantly affected by the proximity of an interface or a nanostructure 
016). Nanostructures made of dielectric materials such as microcavities have been shown to 
induce an increase or a decrease of the fluorescence decay rate Q. Several groups have also 
reported spatial variations of the EM-LDOS in disordered photonic media, based on statistical 
measurements of the fluorescent emission rate from a large collection of dipole emitters (2] 
|8]|9). In such systems, the direct observation of modes strongly localized in sub wavelength 
volumes is crucial as it can provide an important signature of Anderson localization and opens 
the route to cavity quantum electrodynamics in systems which are inherently disordered ifTOl . 
Optical nanoantennas constitute another example of systems in which electromagnetic field can 
be controlled on a nanometer scale, producing an environment which can modify the amount 
of energy emitted by molecules and their direction of emission ifTl ITTI IT2l fT3l IT4l IT31 1X61 . 

A key issue to probe the EM-LDOS experimentally is that the detection process should 
have the same response for all modes. Several methods have been proposed to map the spatial 
variations of the EM-LDOS on photonic nanostructures, among which measuring the thermal 
emission in the near-field flTl . measuring the "forbidden light" signals from the aperture of a 
near-field scanning optical microscope ITT8I . or using a scanning electron beam as a point dipole 
source (19]. The EM-LDOS can be directly inferred from measurements of the spontaneous flu- 
orescence decay rate of a single nanoemitter in its local environment, T = 1/t, where T is the 
fluorescence lifetime of the emitter J2]|6]. An emitter constituted by several randomly oriented 
dipoles probes simultaneously the EM-LDOS in the three spatial directions. One dimensional 
maps of the decay rate were initially obtained by scanning a gold bead near single isolated 
fluorescent molecules fixed on a substrate [20. 21 j. Recently, a decrease of the fluorescence 
lifetime was measured in the reverse situation, when scanning a fluorescent bead across a 250 
nm diameter silver rod, pointing to an increased EM-LDOS due to the existence of plasmonic 
modes on the rod El . 

While a detailed knowledge of the EM-LDOS is clearly required, it is not enough to provide 
a full characterization of a system involving dipoles coupled with plasmonic nanostructures. 
Local changes of fluorescence intensity depend on other parameters such as the radiative and 
non-radiative part of the EM-LDOS and the local field enhancement factor l23l . To characterize 
a plasmonic antenna, one needs at least to measure both the EM-LDOS and the fluorescence 
enhancement factor at the nanometer scale in the near-field of the antenna. 

In this paper, we present the use of a fluorescent 100-nanometers bead grafted at the apex of 
a scanning probe for the near-field investigation of an optical antenna. Two dimensional (2D) 
maps of the decay rate and intensity of the fluorescent bead are obtained with subwavelength 
resolution. The experimental near-field images of fluorescence intensity and EM-LDOS pro- 
vide complementary information that are both needed to characterize the optical response of a 
nanoantenna. A good understanding of these maps is achieved by comparison with numerical 
simulations, which also allows one to analyze the optical response of the fluorescent near-field 
probe used to perform the experiments. 



The decay rate averaged over the orientations of the transition dipole p reads 



r=—\p\ 2 p(r,co), (1) 
3£o« 

where p (r, co) is the EM-LDOS. Hence, measuring the fluorescence lifetime T = 1/T is a direct 
way to probe the EM-LDOS. 

If the molecules are far from saturation, the fluorescence signal reads 



S = C 



^T7(u,0) //i ,„)d£2 



o((D exc )K 2 {(o exc )I inc . (2) 



In this equation, 77 (u, ©//„„) is the apparent quantum yield for a detection in direction u, 
and £2 is the solid angle of the detection objective. The apparent quantum yield is defined 
as rj(u,(Ofi uo ) = rfl(u)/r , with T^(u) the directional radiative decay rate and r = 1/t the 
full decay rate. (0 exc and (Ofi uo are the absorption and emission frequencies of the molecules. 
The constant C is a calibration parameter of the detection (that accounts for transmissivity of 
filters, detector efficiency, etc), a((O ext ) is the absorption cross-section of the bare fluorescent 
beads, 7 (ne is the incident laser intensity and K 2 (co exc ) is the local-intensity enhancement fac- 
tor. In terms of the fluorescence signal, the product F = [/ n Tj(u,0)// ! ,„)d£2] K 2 ((0 exc ) is the 
fluorescence enhancement factor and drives the contrast of the images. 



2. Experimental results 

Active fluorescent probes which allow near-field intensity measurements have been reported 
[ 24 25 26 27 28 29 30 1 , but a scanning probe based on a fluorescent nanobead in which both 
the fluorescence decay rate and the fluorescence intensity are simultaneously measured with 
sub wavelength resolution is still very uncommon J22]. Such a probe constitutes a major advance 
in the field of nano-optics as it provides a direct access to a map of the EM-LDOS and the 
intensity. Its practical realization firstly requires manipulating the bead to graft it at the apex of 
the tip of the scanning probe, and then recording its fluorescence lifetime and intensity at a rate 
compatible with the acquisition of 2D images. The experimental setup which we have designed 
to this aim is sketched in FigQJ. It consists of a homebuilt atomic force microscope (AFM) 
at the tip of which a fluorescent nano-bead is fixed. The AFM is mounted on the stage of an 
inverted confocal fluorescence microscope equipped with a white pulsed laser (Fianum SC450), 
an avalanche photodiode (Micro Photon Devices), and a time correlated single photon counting 
(TCSPC) electronics for fluorescence lifetime measurements (Picoquant HydraHarp 400) 0. 
The scans are performed on the sample stage using a three axes (xyz) piezoelectric assembly 
(Piezosystem Jena). The tip is a pulled silica fiber with an apex of w 100 nm diameter, mounted 
at the extremity of one arm of a quartz tuning fork excited by a dither piezoceramic. Stepper 
motors (Smaract GmbH) allow three dimensional displacements of the tip to perform the coarse 
approach and the fine positioning of the tip within the tightly focused laser spot produced by 
a high numerical aperture (NA=1.4) microscope objective. The tip oscillates laterally in shear 
mode at a height of approximately 10 nm above the surface of the sample thanks to a feedback 
electronics [311 . The optical probe consists of a 100 nm diameter polymer bead filled with 
fluorescent dye molecules (Invitrogen Red Fluospheres) which is grafted at the extremity of 
the silica tip. Fluorescent beads with a diameter of 100 nm are initially spin-coated on a glass 
coverslip with a separation of 5 to 10 jJ.m between them. Wide-field illumination of the coverslip 
is initially used to excite simultaneously the fluorescence of several beads (see FigQJi), which 
are imaged with an EM-CCD camera ( Andor iXon 897) mounted on an extra port of the inverted 
microscope. The silica tip is brought at the glass surface. The grafting of a single fluorescent 
bead on its apex is achieved by moving the bead in contact with the tip while monitoring the 



process in real time with the CCD camera (see movie in supplemental data). FigQJ shows an 
optical fluorescence image of the beads spread on the glass coverslip, one of them being grafted 
on the silica tip visible on the right side of the image. Remarkably, there is no glue involved in 
the grafting which is successful in about 50% of the attempts and can last up to several days. 
Once the grafting procedure is achieved, the nano-bead is brought close to the structure under 
study, the laser (X exc = 560 nm) is tightly focused on the bead in order to excite its fluorescence 
and the sample is scanned. Note that, as the grafting is achieved by lateral displacements of the 
bead towards the tip (see supplemental data), the bead is expected to be laterally displaced with 
respect to the tip apex, as schematically sketched in FigOJ. 

As a first example of the investigation of a nanoantenna, we have studied a linear chain of 
three 150 nm-diameter gold nano-disks separated by 50 nm gaps on a glass substrate 11341 . Ac- 
cording to numerical simulations discussed later in the paper, such a nanostructure is expected 
to produce significant variations of the EM-LDOS and of the fluorescence enhancement factor 
on a scale of several tens of nanometers, which we have observed with the fluorescence scan- 
ning probe. The disk chains are fabricated by electron beam lithography on a glass microscope 
coverslip. Each disk is made of a 2 nm-thick wetting layer of chromium and of a 30 nm-thick 
layer of gold. 

Three signals are simultaneously recorded during the scans as a function of the position 
of the fluorescent probe at the surface of the sample, producing three different images of the 
nanoantenna: (1) the topography (AFM image, see FigfJ^); (2) the fluorescence intensity, which 
corresponds to the integral of the decay rate histogram (Figf2j)); and (3) the spontaneous decay 
rate (Figf2j;), obtained after measuring, at each point of the scan during 1 s, the histogram of 
the arrival time of the fluorescence photons with respect to the excitation laser pulse [2|. This 
long acquisition time is a major difficulty to record 2D maps of the decay rate on a very small 
area with a scanning probe, because the system has to be stable with nanometric precision 
against thermal and mechanical drifts. This explains the low number of pixels (32 x 22) of the 
experimental images. 

The topography of the sample shows three distinct objects with a height of 30 nm each. The 
measured structures seem not to be circularly shaped as expected, but rather elliptical. This 
is because the measured topography is given by the convolution of the topography of the real 
object (three disks) with the shape of the probe used to make the scan, here the AFM tip on 
which a fluorescent bead is grafted sideways. Due to the lateral shift of the bead with respect to 
the tip apex (FigQJ), the topography probe presents an anisotropic shape. The size of the short 
axis in the AFM image of a disk is ruled by the larger object which constitutes the probe, while 
the size of the longer axis is related to the sum of the sizes of the bead and the tip. Hence, the 
topographic image confirms the lateral displacement of the bead with respect to the tip apex 
and allows one to affirm that the size of the tip is of the same order of the size of the bead. 

The geometry of the scanning probe as sketched in FigQJ, is confirmed by looking at the 
fluorescence intensity map shown in FigJJJ}. The contour of the measured topography (Figj2^) 
is reported on this map to guide the eye (dashed lines). The intensity signal decreases in three 
regions, circularly shaped, situated on the upper half of the elliptical contour. Since the fluo- 
rescence signal only comes from the 100 nm diameter dye doped bead, this confirms that the 
three golden disks are scanned twice, once by the bead and then by the tip. The decrease of 
the fluorescence intensity in correspondence of the gold structure is confirmed by the numer- 
ical simulations shown in Figj3] Remarkably, the numerical and experimental intensity maps 
present an almost quantitative agreement, showing both a contrast of about a factor 3 of the 
fluorescence intensity. Note also that the combination of the topography and the fluorescence 
intensity signals allows us to intuit quite precisely what is the position of the bead with respect 
to the tip apex. 
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Fig. 1. a) Sketch of the experimental setup. The active AFM tip is mounted on a piezoelec- 
tric system allowing the positioning of the tip within the laser diffraction limited spot. The 
excitation and fluorescence photons are respectively focused and collected from the same 
high NA objective. The sample can be moved on the XY plane to perform the fluorescence 
intensity, EM-LDOS and topography maps, b) Wide-field fluorescence image of the beads 
spread on the microscope coverslip. The tip can be seen on the right side of the image. One 
bead is grafted on the apex of the tip (see text and movie in supplemental data), c) Artist 
view of the active AFM tip scanning the near-field of the gold nanoantenna. 



Beyond the fluorescence intensity map, a deeper insight of the properties of the electromag- 
netic field at the surface of this nanoantenna is obtained by mapping the decay rate of the 
fluorescent bead, which is proportional to the EM-LDOS (see EqQ]i. Fig|2j; shows a decay rate 
map of the scanned area. It is observed that the EM-LDOS increases by about 30% in three 
regions presenting an extension of about 60 nm each and separated by 100 nm. Two of these 
regions are situated between the gold disks, as confirmed by numerical simulations reported in 
Fig|3}; and discussed below. As in the case of the fluorescence intensity map, numerical and 
experimental data are in almost quantitative agreement regarding the expected change of the 
decay rate in the region between the disks with respect to a region far away from the nanoan- 
tenna. Note that the numerical simulations also show the presence of two other regions with 
an enhanced EM-LDOS, on the external sides of the nanoantenna. Experimentally, we are only 
able to see one of these regions, situated on the right side of the antenna. A possible explanation 
for this is an asymmetry of the gold structure, caused for example by a defect of the lift-off pro- 
cess, that would translate in an asymmetry of the structure of the electromagnetic field on the 
surface of the nanoantenna. Numerical calculations with asymmetric shaped nanoantennas have 
been done and produce similar asymmetries in the EM-LDOS images. However, since the exact 
shape of the nanoantenna is not accessible at the required level of resolution, having an exact 
matching between theory and experiment is a very speculative task and the discussion is there- 
fore limited here to a comparison between the experimental results with numerical simulations 
made on an ideal antenna formed by three regularly spaced circular disks. 




Fig. 2. Experimental results: a) Topography of the sample, b) Fluorescence intensity map. 
c) Decay rate (EM-LDOS) map. The contour of the topographic relief (dashed line), as 
measured by the active AFM asymmetric probe (see text), is reported on the three maps to 
guide the eye. 



3. Numerical simulations: comparison with the experiment 

In order to analyze the experimental results precisely, we have performed exact 3D numerical 
simulations. We use a volume integral method based on the Lippmann-Schwinger equation 

2 

E(r,(»)=Eo(r,©) + ^- J [e(r» - 1] G (r,r',(o)E(r' ! (o)d 3 r / , (3) 

where Eo is the incident field, Go the dyadic Green function of the background (vacuum in 
our simulations) and e(ffl) is the dielectric function of gold l35l . The numerical computation 
is done by a moment method, discretizing the trimer into 5 nm cells (see Figj3^). The Green 
function Go is integrated over the cell to improve convergence ll36l . Computing the electric 
field under the illumination of a source dipole En(r, to) = /io® 2 Go(r,ro, (o)p, we can deduce 



the total Green function G from E(r, co) = /ioft) 2 G(r,r , 0))p. The EM-LDOS is then given by 



^l = ^ImTr[G(ro,ro,c)], (4) 
Po h 

where po is the EM-LDOS in vacuum, given by po = (O 2 / (n 2 c 3 ), ko = (o/c and Tr denotes the 
trace of a tensor. The fluorescence enhancement factor, given by Eqf2] is driven by the product 



T}(u,co f i uo )dD. 



K 2 {(O exc ). (5) 



The apparent quantum yield of the molecules can be deduced from the total Green function 
G ||37l - For sake of simplicity, we do not take into account any intrinsic non-radiative decay 
rate that would account for internal losses. Moreover, in our simulation, we do not integrate 
over the directions u of collection but we only consider the direction orthogonal to the plane 
of the antenna. To compute K(co exc ) we solve Eqf3] under the illumination of a plane-wave 
with normal incidence. We average on two orthogonal polarizations to mimic the fact that the 
incident laser beam in the experiment is not polarized. 

In Figf3] we present the results of our calculations. Every point of each map is averaged over 
100 positions of the emitter randomly chosen inside a 100 nm diameter sphere, to mimic the 
experimental fluorescent beads. The distance between the top of the trimer and the bottom of 
the bead is fixed at d = 20 nm. The excitation and emission wavelengths have respectively been 
set to X exc = 560 nm and Xfi lw = 605 nm. The dimensions of the trimer are the same as the 
experimental ones. 

The general trends observed in the experimental maps are in excellent agreement with the 
simulations. The decay rate exhibits two major hot-spots in the two gaps between the disks, 
and two minor ones on both sides. The fluorescence signal is significantly reduced when the 
molecules are on top of a disk. The contrast of both maps are in nearly quantitative agreement 
with the experimental results, which confirms us in the idea that the order of magnitude of the 
distance between the bead and the trimer is sound. 



4. Numerical study of the resolution of the EM-LDOS map 

One interesting feature of both the experimental and numerical EM-LDOS maps is that both 
seem to exhibit variations on scales well below lOOnm, the size of the fluorescent bead. To 
explain this phenomenon, already observed in 11221 . we compare the contribution to the decay 
rate of the emitters located in the lower half and the upper half of the bead. FigH^-b show 
the EM-LDOS maps averaged over 100 emitters located at different positions inside a lOOnm 
bead for two distances d between the bottom of the bead and the top of the trimer. Fig|4j;-d 
and FigH^-f show the maps averaged over the emitters located respectively in the lower and the 
upper half of the bead. Each map is normalized by the value of the EM-LDOS in vacuum to 
allow for the comparison between the maps. 

A comparison between Fig 0^ and b shows that the smallest details present on the EM-LDOS 
map (such as for example the two EM-LDOS hot spots visible on the right and left sides of the 
nanoantenna) are washed out when the distance of the bead to the sample surface increases. 
This is a known feature in near-field optics. In fact, subwavelength details are evanescent and 
can only be probed by emitters in the very near-field of the system 0. Since these details are 
visible on the experimental map, this study confirms that the real distance between the bottom 
of the probe and the sample surface is of the order of 20 nm. 

Furthermore, a detailed observation of FigJlk, b and c allows us to affirm that the resolution 
of the EM-LDOS map is not limited by the size of the bead. Indeed, the similarity between 




Fig. 3. (a) Top view of the topography of the discretized trimer. Note that the trimer is 
30nm thick, like the one used in the experiment; (b) Numerical fluorescence signal map 
(expressed in arbitrary units); (c) Numerical EM-LDOS map normalized to its value in 
vacuum. 



FigHk and c clearly shows that the EM-LDOS map is driven by the emitters situated on the 
lower half of the bead. The two EM-LDOS hot spots which are visible on the right and on 
the left side of the nanoantenna are smeared out when considering only the contribution of 
the emitters populating the upper part of the sphere. More insight can be given by plotting 
the section of the EM-LDOS maps along the lines drawn in Fig|4^, c, e. The obtained profiles 
are normalized by the maximum value of the corresponding map p max , in order to quantify the 
contrast of each hotspot. They are reported in Figj4g and in Fig|4ji for d = 20 nm and d = 50 nm 
respectively. For d = 20 nm, the lateral hot-spot is clearly resolved when the EM-LDOS signal 
is averaged on the emitters located on the bottom of the sphere or over all the sphere, while it 
is washed out when the signal is averaged over the top of the sphere. Therefore the resolution 
of this detail is clearly due to the bottom emitters. Consequently, the effective resolution is not 
limited by the size of the bead but is smaller and in the case presented in this paper is of the 
order of 50 nm. Interestingly, at d = 50 nm, even if a non-monotonic behavior is observed in 
the profiles reported on Figj4ji), the bottom emitters are too far away from the sample surface 
and the smallest details are washed out. 

5. Conclusion 

In conclusion, the scanning probe described in this paper has demonstrated its ability to map 
both the fluorescence intensity and the EM-LDOS simultaneously with nanometer accuracy. 
While its use was devoted here to the investigation of an optical antenna, it can be used to 
the investigation of the electromagnetic modes on any type of photonic nanostmcture. The 
manipulation of a fluorescent nano-object with nanometer accuracy at the surface of a photonic 
nanostmcture opens interesting new perspectives for studies of quantum electrodynamics, such 
as the investigation of the coupling of quantum emitters with plasmonic (or dielectric) devices, 
the characterization of electromagnetic modes on photonic nanostructures, or the search for 
Anderson localized modes in disordered systems. 
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Fig. 4. Computed normalized EM-LDOS maps for two distances d between the bottom of 
the bead and the top of the trimer. (a-b) Average over 100 emitters randomly located in 
the bead; (c-d) Contribution of the 48 emitters located in the lower half of the bead; (e-f) 
Contribution of the 52 emitters located in the upper half of the bead, (g-h) Section view 
of the maps (a,c,e) and (b,d,f) respectively along the line shown on the maps. Note that in 
this case the EM-LDOS has been normalized by the maximum value of each map p max to 
quantify the contrast of the image. X eX c = 560nm; Ay; i;o = 605 nm. Diameter of the bead: 
lOOnm. 
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